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Abstract. An electron distribution with a temperature 
anisotropy > 1 can lead to the Electron Firehose 

instability (Here || and _L denote directions relative to 
the background magnetic field Bq). Since possible parti- 
cle acceleration mechanisms in solar flares exhibit a pref- 
erence of energizing particles in parallel direction, such 
an anisotropy is expected during the impulsive phase of a 
flare. The properties of the excited waves and the thresh- 
olds for instability are investigated by using linearized 
kinetic theory. These thresholds were connected to the 
pre-flare plasma parameters by assuming an acceleration 
model acting exclusively in parallel direction. For usually 
assumed pre-flare plasma conditions the electrons become 
unstable during the acceleration process and lefthand cir- 
cularly polarized waves with frequencies of about ^ \ flp\ 
are excited at parallel propagation. Indications have been 
found, that the largest growth rates occur at oblique prop- 
agation and the according frequencies lie well above the 
proton gyrofrequency. 
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1. Introduction 

Particle acceleration is a phenomenon occurring at many 
different sites throughout the universe. An important ex- 
ample of particle acceleration are solar flares, offering a 
wide range of observations that allow one to probe elec- 
tron and ion acceleration. It is now widely accepted that 
the hard X-ray emission observed by various spacecraft 
reflects the energization of almost all electrons in the flar- 
ing plasma to energies up to ^ 25keV. These observations 
and the observed magnetic fields encompassing the solar 
fiare suggest that most of the dissipated energy is released 
by restructuring the magnetic field, e.g. magnetic recon- 
nection events. 

During the impulsive phase of the flare, when the most 
powerful energization takes place, electrons must be ac- 



celerated to mean energies of about ^ 25keV at a rate of 
about 10^^ electrons per second in order to sustain the 
observed intensity of the hard X-ray bursts. Taking the 
impulsive phase of a flare to last about 10s and assum- 
ing a n electron density of abo ut 10^°cm~^ (Moore &Fung 
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I972| ; Vaiana & Rosner |1978| ) , the bulk energization must 
process a coronal volume of at least lO^^cm^. Thus ener- 
gization must affect a large fraction of the electron popu- 
lation in the flaring region. 

In view of this background we want to briefly describe 
the processes which may be responsible for particle accel- 
eration in solar flares. For a detailed review of possible ac- 
celeration processes in impulsive solar flares see e.g. Miller 
et al. ( p^ ) and CargiU (|l99|). 

1) Shock Acceleration: There are two types of shock ac- 
celeration. The one referred to as 'shock drift acceleration' 
involves the shock electric field that refiects and acceler- 
ates the particles moving along the shock surface. Since 
this mechanism is only effective when the shock normal 
approaches an angle of 90° to the background magnetic 
field, either the gained energy or the particle flux is very 
limited. It seems to be unlikely that this mechanism is 
responsible for the large number of accelerated particles 
in solar flares. The second kind of shock acceleration is 
called 'diffusive shock acceleration'. In this process the 
particles cross the shock-front several times, interacting 
with scattering centers on both sides of the shock. In the 
rest frame of the shock these centers approach each other 
and the particles systematically gain energy. This kind of 
acceleration process requires a certain initial velocity in 
order to become effective. The ion velocity has to exceed 
the Alfven speed Va = Bq/ y^^irp while the electrons must 
have velocities at least above ^JrrvJm^K. This has been 
called the 'injection problem'. 

2) Acceleration by parallel electric fields: Direct accel- 
eration by electric fields depend on its strength compared 
to the Dreicer field Ed = (27re3nc In A)/(fcBT). If > E^d 
most electrons and ions gain energy. If iJ < E-a only elec- 
trons in the high energy tail of the velocity distribution 
function will be accelerated. The limitation in both cases 
is the maintenance of overall neutrality of charge and pre- 
existing current in the acceleration region. 
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3) MHD turbulence: This acceleration mechanism oc- 
curs when particles interact many times with randomly 
moving MHD waves. Due to a slight overplus in head-on 
collisions the interaction results in an energy gain for the 
particle. As in the shock acceleration model, the acceler- 
ation by MHD turbulence suffers from an 'injection prob- 
lem': Thermal ions and electrons cannot resonate with 
MHD waves for typical solar pre- flare conditions. 

A solution for this problem is the assumption of MHD 
turbulent cascades that channel the energy residing in 
the MHD turbulence to smaller scales and into the re- 
gion where interaction with thermal particles is possible. 
A realization of th is sce nario is prop osed i n Miller ( 1991 ), 
Miller & Roberts (|l995D and Miller (|l997|) . An MHD tur- 
bulent cascade transfers the energy from large scale MHD 
waves to smaller scales where the energy may be absorbed 
by the particles. The mechanism that dissipates the wave 
energ y int o the particles is transit-time damping (Fisk 
197^ ; Stix 119921) . It is is basically a resonant Fermi accel- 



eration of second order. Only particles in resonance with 
a low-amplitude MHD wave are affected. The resonance 
condition is the usual / = (or Landau) resonance given 
by a; — « 0. As the particles can only gain energy 
in the direction parallel to the background magnetic field, 
the temperature in parallel direction increases. 

A preference for acceleration along the background 
magnetic field is a common feature of the acceleration 
models mentioned above. The velocity distribution thus 
becomes more and more anisotropic during acceleration. 
If energization in parallel direction is from a thermal level 
of some 0.1 keV to 20 keV or more but the perpendicular 
temperature remains constant, the anisotropy is substan- 
tial. The free energy residing in parallel direction may give 
rise to growth of plasma waves. 

For Ti? > and high beta plasmas, HoUweg & Volk 



( |1970D and Pilipp & Volk ( |1971| ) have proposed the Elec- 
tron Firehose instability. This instability is an extension to 
higher frequencies of the (MHD) Firehose instability, orig- 
inally mentioned by Parker ( |1958| ) . While the Firehose in- 
stability is of a completely non-resonant nature, the Elec- 
tron Firehose instability involves non-resonant electrons 
but resonant protons. 

For large anisotropy of the electron distribution, the 
electrons become also resonant. This instability is de- 



scribed in Pilipp & Benz (1977) and is called the Resonant 
Electron Firehose instability. 

Having been applied to a variety of problems, the Elec- 
tron Firehose instability has not been considered to occur 
during electron acceleration in solar flares. 

Here we investigate the threshold for growth of plasma 
modes resulting from acceleration and infer a prediction 
for the evolution of the distribution function in velocity 
space with respect to conditions expected to occur in solar 
flares. 

We assume that no significant instability of Langmuir 
waves occurs. This is suggested by the following argument: 



If a large fraction of the available energy normally did 
go into Langmuir waves, we would expect always a radio 
signature orders of magnitude higher than ever observed 



during the impulsive phase (Benz & Smith 1987). 

Section 2 describes the techniques used to solve the 
dispersion equations. In the following section the results 
obtained are shown and the thresholds of the instability 
are presented. In section 4 we discuss the effect on the 
acceleration of electrons and conclude this work. 

2. Method 

Consider electromagnetic transverse waves of the form 
exp {ikx — iut) propagating in the direction of the back- 
ground magnetic field in an electron-proton plasma. The 
plasma dispersion equation can then be written as 



det ( {Ik^ - kk)— - e(w, fc) ) = 0, 



(1) 



where, according to linearized kinetic theory, the dielectric 
tensor e{Lj,k) is given by 



e(w, fc) = 



J 71 — — OO 




(2) 



The gyrofrequency of the jth species is given by ilj = 
qjB/{cmj) and Wp denotes the plasma frequency defined 

as t^p = y^Ej^ = y^Ej ^Trrij-gl/TOj. The tensor J] is 
given by the matrix 



n 



(3) 



■ nilj J jf 
iv±Vl\JnJn 



'iv±V\lJnJ'^ 

(w||J„)' / 



where the argument of the Bessel function J„ is k±v±/Qj. 
/j (w|| ,v±) in equation (||) denotes the zero order distribu- 
tion function in velocity space of the particle species j. In 
order to obtain full solutions of this equation, the com- 
puter code WHAMP (Ronnmark 1982| ) has been applied 
to the problem. The usage of this code has been facili- 
tated by programming an interface for the programming 
language IDL. It is called IDLWhamp and provides the 
user with a comfortable tool to input parameters to the 
code and for management and visualization of the results. 

According to the capability of the WHAMP code, the 
most general form of the particle distribution function 
fji^\\,v±) for each species j is given by 
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2TTV 



■ cxp 



A, 



3th 



■ exp 



exp 




(4) 



2a2jV 




This is the original notation used in Ronnmark (1982) 
beside the choice of the thermal speed to be Wj^^ = 
[k-QTj /m.j). The aij parameter is the anisotropy aij — 
T]_/T^t^ of the j-th distribution function. and a2j define 
the depth and size of a possible loss-cone. 

We assume that the electron velocity distribution func- 
tion can be described by a bi-maxwellian with different 
temperatures in parallel and perpendicular direction with 
respect to the background magnetic field. Hence for our 
problem the parameters and a2j were set to unity in 
equation (^. 

Taking into account the uncertainties in the acceler- 
ation region including a possible pre-heating mechanism 
altering the pre-flare plasma conditions, we do not want 
to restrict our work to the parameters of a particular sce- 
nario. According to Pallavicini et al. ( 1977] ) reasonable 
ranges in the acceleration region of an impulsive solar flare 
would be ~ 100 — 500G for the background magnetic field, 
« 10^-10^^cm~^ for the number density and w 10^-10''K 
for the temperature of electrons and protons. 




CO. 



Fig 

sen parameters are T' 



1. A typical plot of the dispersion relation. The cho- 
Tl = TP = lO^K, ^ 



20,71, 



5 • 10^cm~'^,Bo = lOOG. The real part of the frequency 
Wr and the growth rate 7 are normalized to the proton gy- 
rofrequency |f2p|. The parallel wave vector is normalized 
to the proton inertial length. The whole branch is lefthand 
circularly polarized. 



3. Results 

3.1. Electron Firehose Instability 

The only mode exhibiting significant growth rates in our 
calculations is a lefthand circularly polarized wave which 
was identified to be the Electron Firehose instability. This 
mode evolves out of a stable righthand polarized whistler 
wave at small anisotropy. With increasing anisotropy, the 
frequency Wr is shifted so that, with the convention uj-c > 0, 
the mode becomes lefthand circularly polarized at A;||Bo in 
the unstable regime, cf. section 7 in Gary ( 1993| ). A typical 
dispersion relation is plotted in Fig. ^ By introducing the 
resonance factor 



^\h\v. 



(5) 



the values for the protons and electrons arc found to be 
\C,~\ ~ 1 and 1, demonstrating resonance for the 

protons and non-resonance for the electrons. 



According to HoUweg & Volk (1970) there are also 



right hand circularly polarized modes, which can become 
unstable for this extension of the Firehose instability. 
These modes have been found, but the growth rates are 
smaller than the growth rates of the lefthand polarized 
modes described above. 
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Fig. 2. Resonance factor |Cp | of the protons versus the 
electron anisotropy ^ for the fastest growing modes. The 
larger |Cp^|, the smaller is the fraction of protons in reso- 



nance. 



As the instability first appears, the phase velocities of 
the resonant waves are near the peak of the proton dis- 
tribution. Fig. H displays a plot of the proton resonance 
factor (H) for the fastest growing modes versus electron 
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and the usual parallel plasma beta by 



anisotropy. With increasing anisotropy less protons are 
resonant and the resonance factor increases. The change 
in the fraction of resonant protons is mirrored in the ex- 
cited frequency range. As depicted in Fig. || the unsta- 
ble frequency range grows to a maximum value at an 
anisotropy of about ~ 12, coinciding with the max- 
imum value of the resonance factor at |Cp"| ~ 0.57 of 
the protons (cf. Fig. ^). As the resonance factor decreases 
again, the excited frequency range becomes narrow around 



(7) 



where the connection between these two quantities is given 
by the temperature anisotropy T^/Tj^. 

According to HoUweg & Volk (197C) the instability 



criterion for the Electron Firehose instability may be ap- 
proximated by 



in. 



1. This narrow range is in itself evidence for the ^ P\\^c < 0, 



(8) 



resonant character of the instability. 



3.2. Instability Threshold 

In this section we present the calculated threshold for lin- 
ear growth of L-mode waves excited by the Electron Fire- 
hose instability. 

The initial plasma is assumed to be maxwellian with 
temperatures Tj_q — T||g — Tq perpendicular and parallel 
to the background magnetic field for both plasma species, 
the electrons and the protons. Taking into account an 
acceleration mechanism for the electrons acting only in 
parallel direction, the perpendicular temperature remains 
constant throughout the whole acceleration process, i.e. 

= T^g. In order to investigate the condition of the pre- 
flare plasma for the Electron Firehose instability to occur 
during the acceleration process, the initial plasma parame- 
ters have to be connected to the actual plasma parameters 
during the acceleration. With the assumptions above, this 
can be done by defining an initial parallel plasma beta, 
/3|jg, via the perpendicular plasma beta 



(6) 



where the anisotropy factor is defined by Ac = 1 — T1/Tp 
As one can see from inequality (^, the instability 
threshold does not depend directly on the parameters 
no, Ti?, Bq, but only on the resulting This indepen- 



dence is also reproduced with the numerically obtained 
data. For our purpose, the plasma is therefore fully de- 
scribed by the plasma beta. 

In Fig. I the function 7max(T||=/T|) is plotted for five 
different values of /Jjjg. The maximum growth rate of the 
instability steeply raises at the threshold of the instabil- 
ity and flattens for larger anisotropics, where 7max/|^^p| 
approaches unity. 

From these results, the contour of zero growth rate in 
the Ac — plane has been derived (cf. Fig. ||). The dis- 
crepancy between the analytically derived relation (^) and 
the numerically obtained values is due to the approxima- 
tion used in HoUweg & Volk (|1970| ). 

According to inequality (^, instability cannot occur 
for a parallel beta smaller than unity. Due to the devia- 
tions of the approximation mentioned above, this limit is 
shifted to a value of ~ 1.6 (cf. Fig. |^). 

In order to investigate the necessary properties of the 
pre-flare plasma for the Electron Firehose instability to 
occur, it is the initial plasma beta that is of interest. Fig. || 
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Fig. 5. Threshold for the Electron Firehose instability in 
anisotropy factor vs. parallel electron beta. The scaling 
of both axes is logarithmic. The dashed curve shows the 
instability limit according to equation (^) . The dotted line 
represents a fit to the numerically obtained values which 
are depicted by diamonds. The areas above the respective 
lines are the unstable regions. 
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Fig. 6. The same plot as Fig. ^ but this time Ac has been 
plotted vs. the initial electron beta. Again, as in Fig. ||, 
the dashed curve shows the instability limit according to 
equation and the dotted line represents a fit to the 
numerically obtained values. 



depicts the same plot as Fig. ^ but this time the anisotropy 
factor Ac has been plotted versus the initial plasma beta 
The dotted line in both figures represents a fit to 
the numerically obtained values and is an extrapolation 
to a broader range of beta values. The negative at the 
Ac 1 limit is an artifact of this extrapolation. 

The values of the initial plasma beta for the Elec- 
tron Firehose instability to occur at considerable values of 
Tji'/T^ are well within the range of usually assumed pre- 
flare plasma parameters. For example, an initial plasma 



beta of « 0.05 can be realized by assuming pre-flare 
plasma parameters of rie = 5 • 10^°cm~^, Tq = 3 • 10^ K 
and Bq = lOOG. This plasma becomes unstable at an 
anisotropy of T^^/Tj_ « 32. 

According to the acceleration model via transit-time 
damping, this is a reasonable value for the anisotropy to 
occur during the acceleration process (Lenters & Miller 
1993). 



3.3. Influence of Anisotropic Protons 

If we assume the protons to be heated by the same or 
a similar mechanism, it is to be expected that they will 
grow anisotropic in the same way the electrons do. Hence, 
we also have investigated the influence of anisotropic pro- 
tons and briefly discuss the effect of an additional proton 
anisotropy on the instability. 

Consider a plasma with anisotropic electrons and 
isotropic protons that is already unstable to the Electron 
Firehose instability. When the protons are anisotropized 
by increasing the temperature in parallel direction, more 
and more become resonant to the L- waves, non-resonantly 
excited by the electrons. As shown by Hollweg & Volk 
( 1970| ), the protons are damping these waves. Hence, it 
is to be expected that the resulting growth rate of the L- 
waves decreases as the proton anisotropy is increased. This 
expectation has been verified by numerical calculation. 

Moreover, the protons are heated by absorption of the 
excited waves at the expense of the electrons (Pilipp & 
Volk |l97l| ). According to Kennel & Petschek ( [1966D this 
scatters the protons to higher perpendicular velocities and 
hence, destroys or inverts the parallel proton anisotropy. 
It inhibits the growth of the electron anisotropy and may 
complicate the acceleration to higher energies, but in- 
creases the bulk energy of the protons. The energy transfer 
from the electrons to the protons via the Electron Fire- 
hose instability could be responsible for the proton ener- 
gization, which is a problem in the transit-time damping 
scenario (Miller |1998| ). 

If the protons are anisotropic, there is an additional 
righthand polarized wave mode. This mode is the ex- 
tension of the lefthand Electron Firehose mode toneg- 
ative frequencies. According to Hollweg & Volk ( 1970 ) 
this mode has real frequencies mainly below the proton 
gyrofrequency and is in resonance with the protons. As 
the anisotropy of the electrons increases, this righthand 
polarized mode becomes less and less significant. 

The proton cooling through the instability of the right- 
hand mode competes with the heating by the lefthand 
mode and it is not yet clear what the net energization of 
the protons will be. 



3.4. Oblique Propagation: Preliminary Results 

Sample calculations in oblique directions indicate an un- 
stable branch of solutions that grows faster than the par- 
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Fig. 7. Frequency and the according growth rates vs. the 
propagation angle with respect to the background mag- 
netic field. The dashed lines indicate the branch of the 
parallel Electron Firehose instability. The solid lines show 
the oblique mode that exhibit much faster growth than 
the Electron Firehose instability. The small plot is an en- 
largement of the region, where the crossing of the growth 
rates occurs. The plasma parameters are the same as in 
Fig. 1 

allel Electron Firehose mode. This mode is stable at par- 
allel propagation and is also lefthand circularly polarized. 
Fig. ^ depicts frequency and maximum growth rate versus 
the angle 6 of both modes. The dashed lines represent the 
branch, that is excited by the Electron Firehose instability 
at parallel propagation. At a propagation angle of about 
10° with respect to the background magnetic field, the 
growth rate of the oblique mode becomes larger than the 
growth rate of the parallel mode. Hence, the determining 
mode for instability thresholds is the oblique mode rather 
than the parallel Electron Firehose mode. 

As calculations have shown, not only the growth rate 
of the oblique mode is larger than in the parallel case, but 
also the instability threshold may be lower with respect 
to anisotropy. Plasmas being stable with respect to the 
parallel mode exhibited instability to the oblique mode. 
Therefore, the thresholds for instability derived in the sec- 
tion above can be considered as upper boundaries. 

Lefthand circularly polarized oblique modes seem to 
have never been considered in connection with the Elec- 
tron Firehose instability. The further investigation of these 
modes is the subject of ongoing work. 

4. Discussion and Conclusion 

Numerical solutions of the dispersion equation for lefthand 
circularly polarized electromagnetic waves, propagating 
parallel to the background magnetic field, have shown that 
the Electron Firehose instability, usually considered as a 
'high-beta plasma' instability, must be expected in coro- 



nal plasmas which are processed by an acceleration mech- 
anism with a preference in parallel direction. The distri- 
bution function of the electrons in velocity space has been 
represented by a bi-maxwellian with temperatures and 
T|, perpendicular and parallel with respect to the back- 
ground magnetic field. The protons have been assumed to 
be isotropic and in thermal equilibrium. 

Considering the uncertainty in the pre-flare condi- 
tions, we have investigated instabilities in a broad range 
of plasma temperature T, density n and background mag- 
netic field Bq. For these plasmas, it was found that the 
Electron Firehose instability occurs at anisotropies that 
must be expected for an acceleration mechanism acting 
predominantly in parallel direction and being capable of 
producing the observed electron bulk energization. 

The unstable parallel modes that have been found are 
lefthand circularly polarized and non-resonantly excited 
by the electrons, but partially cyclotron resonant with the 
protons, which absorb the wave energy. Hence, energy is 
transferred from the electrons to the protons. Assuming 
the density and the magnetic field to be constant dur- 
ing the acceleration process, there is a limiting electron 
temperature in parallel direction that cannot be exceeded 
without loosing energy to the protons via the Electron 
Firehose instability. The Electron Firehose instability may 
thus inhibit the electron acceleration process and limit the 
reachable energies. 

At angles 9 7^ with respect to the background mag- 
netic field, an oblique mode has been found, that exhibits 
even larger growth rate than the parallel Electron Firehose 
instability. This mode is also lefthand circularly polarized. 
The properties of the oblique mode open new aspects on 
the Electron Firehose instability and its thresholds. 

Taking anisotropic protons into account, the lefthand 
mode may extend to negative frequencies. They corre- 
spond to the righthand circularly polarized mode reso- 
nantly excited by the protons. Due to the cooling effect of 
this instability it is not yet clear to what the net energiza- 
tion of the electrons and the protons will amount. This 
topic will be the subject of future particle simulations. 
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